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Abstract. 

The analysis which assumes that tick by tick data is linear may lead to wrong 
conclusions if the underlying process is multiplicative. We compare data analysis done 
with the return and stock differences and we study the limits within the two approaches 
are equivalent. Some illustrative examples concerning these two approaches are given. 
Actual data is taken from S&P 500 stock cash index. 



1. Introduction 

One of the most important problems in mathematical finance is to know the probability 
distribution of speculative prices. The first approach to the problem was given by 
Bachelier in 1900 when he modelled price dynamics as an ordinary random walk where 
prices can go up and down due to a variety of many independent random causes. 
Therefore, the distribution of prices has to be Gaussian |TJ due to the Central Limit 
Theorem: the sum of independent, or weakly dependent, random disturbances, all of 
them with finite variance, results in a Gaussian random variable. 

Despite Bachelier's very early interest in stochastic modelling of stock prices, 
research on this topic is not again noticeable until 1930 's. A renewed regard on 
financial markets appeared in the embryo school of American economists highly skilled 
in mathematics and statistics. In an ideal and theoretical framework, they believed 
that market was perfect in the sense that one cannot forecast future price changes 
based on past history alone. Therefore, they conclude that price changes have to be 
uncorrelated, and follow a Gaussian random process thus obeying the Central Limit 
Theorem enunciated above. 

At that time, the main research in finance was addressed to test those theoretical 
hypothesis on real markets. In 1953, Kendall analyzed several American markets 
observing there that price changes behave like wandering series and discovering 
correlations in the price movements time series ||. All of this partially contradicted 
the economic theory since, although Kendall confirmed the random nature of stock 
evolution, he also found correlations which were unacceptable in an ideal and perfect 
market framework. 
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Few years later, Osborne tried to fit data with a Brownian motion model and looked 
for the form of the price changes empirical distribution ||. Price movements modelled 
as a random walk implied that price can be negative with non-zero probability. In order 
to avoid the complications posed by the fact that stock prices must have a lower bound, 
he proposed to take the logarithm changes of prices instead of the price changes, i.e., 
\n[S(t + A)/ S(t)] instead of S(t + A) — S(t) . In this case, there were no need of limiting 
the process into a positive region. Osborne confronted new model with real markets 
and observed that the new variable (call it stock return, R(t) = \n[S(t)/ So}) evolved 
as a Gaussian random walk and had its increments uncorrelated. This feature was in 
concordance with the perfect market hypothesis, and gives more importance to the way 
data is taken from markets. 

This historical introduction do not want to review the existent market models nor 
present a new model. From sixties to nowadays, we have much larger time series recorded 
from markets, and computers have allowed us to register the whole stock movements, 
the so-called tick by tick data or high-frequency data. People concern in financial 
markets have dramatically increased and there is a strong demand of high precision in 
the description of the speculative prices dynamics [[| . Data analysis in financial markets 
has thus become a relevant issue, and taking good estimators for checking economic 
theory and market models is an essential but delicate task. Kendall's and Osborne's 
works exemplify how important is to think over which data shall we handle and how we 
manipulate this data. Indeed, the way data is manipulated may lead to diverse and, in 
critical cases, contradictory conclusions. 

More precisely, the purpose of this paper is to consider the way we analyze the 
financial market data. We focus in the differences between taking stock price and 
return increments, and we show their accuracy and range of validity for estimating 
parameters describing the market. Our intention is to display risk of obtaining wrong 
conclusions when we operate in an inadequate framework with the historical time series. 
For instance, we see the different aspect adopted by the probability distribution in 
Figures [I] and || where we respectively take stock differences or return differences time 
as a data source. Database to make the comparison is the Standard & Poor's 500 stock 
cash index ranging from 1983 to 1999, and we fulfill the assertions with a simple market 
model, the multiplicative Gaussian model which posseses properties also assumed in 
more sophisticated and realistic market models. 

The paper is divided into six sections. In section 2 and 3 we present several functions 
related to the return and stock differences, respectively. We study their properties in 
general but also for the Wiener process particular case. In section 4 we show and 
compare the empirical probability distributions for S&P 500 cash index for the stock 
and return differences. Section 5 concentrates on the estimators of the first and second 
moments. Conclusions are drawn in section 6. 
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Characteristic functions are very useful for evaluating moments of any order. We depart 
from the stock return defined as the logarithm of the stock price, i. e., R(t) = \n[S(t) / S ] , 
and we will thus calculate expressions related to the return stochastic variable. We will 
do it in order to obtain these expressions in terms of the return characteristic function ||. 

The characteristic function is derived from the return probability density function 
(pdf) and defined as follows 

/oo 
drj" r p R (r,t\0), (1) 
-oo 

where the conditional density is defined as pn(r, t\0) = p R (r,t\r — 0, t — 0). 

Return differences and the particular case assuming return evolution to be driven 
by a Wiener process are also studied. All these calculations are done with the demand 
that process is Markovian and homogeneous. 

2.1. Some expressions for the return stochastic variable 

Let us present some results related to the return variable using equation ([I]). For 
instance, the first moment of the return is 

(R(t)) = -id UJ Mu,t\0) (2) 

and the second moment is 

{R{tf) = -dij R {u,m\ (3) 

We can also derive the variance 

Var[i2(t)] = (R(t) 2 ) - (R(t)) 2 = - 0*0", *|0) + [d u M«>, W n , (4) 

where equations (0) and (|3|) have been taking into account. 

On the other hand, we can obtain the joint probability density function at two 
different times with the condition that underlying process for the return is Markovian, 
i.e., 

PR(r 1 ,ti 1 r 2 ,t 2 \0) = PR(ri,ti\r2,t 2 )pR(r2,t 2 \0), (5) 

whenever ti > t 2 . But, if we also impose that process is homogeneous in time and 
return, i.e., 

PR(ri,ti\r 2 ,t 2 ) = p R (n ~ r 2 ,h - t 2 \0), (6) 
We can go further and see that 

/OO /"OO 
dn e iwiri / dr 2 e iw2r2 p R {n, t ± - r 2 , t 2 \0) 
-oo J—oo 

/oo roo 
dn / dr 2 ^ M r 2) Pr{ti _ f2)ti _ t2 | 0) pj? ( r2;t2 | ), 
-oo J—oo 

and we thus obtain the joint characteristic expressed as a product of two characteristic 
functions. That is: 

(f>R(wi,tx;w 2 ,t 2 \0) = <£r(o>i,£i -t 2 \0) (f) R (uJi +u 2 ,t 2 \0). (7) 
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Indeed, the correlation function for the return can be written in terms of the joint 
characteristic function. Using equation (|7|), we can derive the correlation function which 

reads 



(RfaRfo)) = -d 2 6 R (^h;u; 2 ,t 2 \0) 



Wl ,UJ2=0 

d,„6n(uj, U — U\0) 

oj=0 



Taking into account equations (0) and (|3]), we write an expression for the correlation 
function in terms of the first and second moments of R. Thus, 

(R(tx)R(h)) = (R(t 2 ) 2 } + (i^(t 2 ))(i^(t 1 - t 2 )) {t x > t 2 ). (8) 

It is usually defined a coefficient which evaluates the degree of correlation between 
a pair stochastic quantities J7j. The coefficient p here defined is enclosed between the 
interval p = [—1,1]. In case that p — 0, it is said that the pair of stochastic quantities 
are uncorrelated. And in any other case, we shall say that quantities are correlated 
and that the correlation is positive or negative according as p > or p < 0. When the 
coefficient raises one of its extreme values, it is said that one quantity is a linear function 
of the other, and the two quantities vary in the same linear sense, p = 1, or in inverse 
sense, p — — 1. For the case referred to the stochastic return variable, the correlation 
coefficient reads 

(*(«,)*(«,)) - WW) (9) 

IVar[R(tx)} V&r[R(t 2 )} 



We can simplify this expression with the help of equations (|) and @ 
( , , ^ _ W2) 2 ) + W 2 )){R{h-t 2 )) - {R{t 1 )){R{t 2 )) 

y/Va^R^] Var[ J R(t 2 )] 
We need to know an equivalent expression for (R(ti — t 2 )) 

(R(ti-t 2 )}= / da a pr{q, t x - t 2 |0) 



—00 
00 



/oo roo 
dr 2 / da a p R {a,t x - t 2 |0) p R (r 2 ,t 2 \0), 
-00 J — 00 



where we have only added an expression which value is one due to the fact p(r 2 ,t 2 \0) 
is normalized. Taking into account that process is Markovian and homogeneous whose 
definitions are given by equations ©-(H) and doing the change of variables a = r\ — r 2 , 
we have 

/oo roo 
dr 2 I dr x (n - r 2 ) Pfl(ri,ti|r 2 ,t 2 ) Pr^M^) 
-00 J — 00 

roo 

dr 2 \ dr x (n - r 2 ) p R {ri, t x ; r 2 , t 2 \0). 



—00 



Hence, 



(Rih - t 2 )} = (Rfa) - R(t 2 )} = (Rih)) - (R(t 2 )}. (10) 
After simple manipulations we finally obtain 

o(t t)~ I Var ^^)] (U) 

p{tl ' t2) -\v ai [R(t 1 )Y 
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2.2. The stock return difference 

However, the variable in which we are specially interested is the one referred to the 
return differences. We define a new stochastic variable called stock return differences 
by 

W(t;r) = R(t + r) - R(t). (12) 

This stochastic variable has the same pdf and, therefore, same moments and 
correlation function as the return. Let us show this. From equation (0), we see that 

/•oo roc 

Pw(w,t;r)= / dr dr' 5[w — (r — r')] p R (r,t + r;r' \t\0) 



— oo J — oo 

OO POO 



oo J — oo 



dr / dr' 8[w — (r — r')] p R (r — r', r|0) PR(r', t\0), 



where we have taken into account equations (|5|) and (|^). Implementing the delta function 
in the inner integral of this expression we get 



p w (w, t; r) = p R (w, r|0) / dr p R (r - w, t\0), 

J—oo 

and since pFt(r, t\0) is normalized we finally obtain 

p w (w, t; t) = p R (w, r|0), (13) 

which shows that pw is identical to the return pdf p R . Observe that distribution is only 
function of the time difference r and does not depend on time t. Therefore, we can 
take expressions presented in equations (PD (|TT]) and replace t by r in order to give the 
equivalent expressions for the return differences W(t; r). 

We now study the autocorrelation between the variable W evaluated at distinct 
times t and t' > t + r. Thus, 

[Wit] r)W(t'; r)> = (R(t + r)R(f + r)> + (R(t)R(f)) 

- (R(t + T)R(t')) - (R(t)R(t' + r)>, 

where we have decomposed the function W in terms of the return R. Equation (|8]) gives 
us the value of each autocorrelation. Taking into account the requeriment that t' > t+r 
and after simple manipulations, we get 

(W(t; r)W{t>- r)> = (R(t + r))(R(t' - t)) + (R(t))(R(f - t)) 

- (R{t + r))(R(t' -t-r))- (R{t))(R(t' + r- t)). 



We can go one step further with the help of equation (10) which gives the the time 
invariant property for the first moment of the return. In terms of the return differences, 
the autocorrelation reads 

(W(t-r)W(t';T)) = (W(t-r))(W(t';r)), 

and thus see that the W(t; r) quantities are uncorrelated. 

In addition, we observe that the returns increments and the return itself are 
uncorrelated stochastic variables since the correlation function is 

(W(t;r)R(t)) = (R(t + r)R(t)} - (R(t) 2 } = (R(r)) (R(t)) , 
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where equation (|8|) is been applied. And we may write 

{W(t;r)R(t)) = (W(t;r))(R(t)). (14) 

thus being zero their correlation coefficient defined in equation^). We therefore conclude 
that W(t; t) and R(t) are uncorrelated. This is not surprising since Markovian condition 
and homogeneity leads directly to the statement that W(t; r) and R(t) are independent 
stochastic variables. Indeed, if we take the joint pdf, we can first implement Markovian 
condition (||) 

p R (r', t'; r, t) = p R (r' -r,t' - t)p R (r, t), 

and, afterwards, homegeneity @ let us write joint pdf in the form 

p R {r', t'; r, t) = p w {w, r)p R (r, t). 

From last equation, we there see that joint pdf becomes the independent product of two 
other pdfs which are, in effect, the distributions of W(t, r) and R(t). 



2.3. A simple model for the return: The Wiener process 



Let us now study these expressions for the Wiener process with drift /i and diffusion 



coefficient a, a well-known market model in the literature 



PR{r,t\0) 



■ exp 



(r - fxtf 
2aH 



The conditional pdf reads 

(15) 



V2Tia 2 t 

and the characteristic function defined in equation ([]]) is 

4> R (w,t\0) = exp (ifiu - ^a 2 u 2 ^j t . (16) 
From this we see that 

(R(t)) = pt, (R(t) 2 ) = a 2 t + /j 2 t 2 , Vax[R(t)] = o 2 t. (17) 

Moreover, we can also calculate the joint characteristic function obtained with the help 
of equations (0) and (|16D. Thus, 

w 1 ,ti;w 2 ,t 2 \0) = exp j(V^! - T^^i) (*i _t 2)| 



x exp 
exp 



-a 2 (uJi + uj 2 ) 2 



1 



iyu(co>iti + u 2 t 2 ) - ^a 2 (oo 2 t l + >J\t 2 + 2uiu 2 t 2 ) 



(18) 
(19) 



Then, the correlation function according to equation ([8]) is 

{R{t x )R{t 2 )) = a 2 t 2 + fiH 2 tu 

and the correlation coefficient can be obtained with equation ( |TT] ) once we know the 
variance of equation (li~7|) . Hence, 

h 



p{h,t 2 



(20) 
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where since t% > t 2 the coefficient p is positive and less than 1. 

As shown in general, the expressions for the return differences in the Wiener case 
are equivalent to those presented for the return. We note that, for this case, return 
differences and return stochastic variables are also uncorrelated (see equation ([14])). 



3. The stock share price 

The purpose of this section is to derive the functions related to the stock price stochastic 
variable S(t) in terms of the results obtained for the return R(t). The section is 
analogous to the one of the return variable but we now implement the same equations 
to the stock price stochastic variable. We will also study price differences stochastic 
variable and the particular Wiener process as a market model. 



3.1. Some expressions for the stock price stochastic variable 
We remember the relation between the return and the stock price: 
S = S e R . 

Our intention is to derive several expressions related to the stock price and write them, 
at the end, in terms of the characteristic function of the return. Therefore, the first 
moment and the second moment are respectively 

(5(t)) = 5 (e i?w )=5oM-i,t|0), (21) 

and 

(S(t) 2 ) = S 2 (e 2 ^)) = S 2 <M-2i,t|0). (22) 

The variance thus reads 

Var[S(t)] = S 2 [(e 2R ^) - (e R ^) 2 ] = S 2 [fo(-2i,t|0) - R (-i,t|O) 2 ] . (23) 

In addition, the characteristic function for the stock price can also be obtained as 
a sum of characteristic functions for the return. That is: 

roc roc 

<Ps(w,t\S ) = ds e^ s p s (s,t\S ) = dr j" Soe p R (r,t\0) 

JO J-oo 

= ± ^T 1 f dr P*(r, * 10) e- = £ ^ M-™, f |0). (24) 
^o n - J ~°° ^o nl 

We also see from equation(^) that, when t\ > t 2 , 

/oo roo 
dr x / dr 2 S 2 e r ^ ^(r 1; t i; r 2 , t 2 |0) 
-oo J —oo 

= S 2 M-hti-t 2 \0) <M-2i,t 2 |0), (25) 

from where we can write autocorrelation in terms of first and second moments if we take 
into account equations ( pl|) and (|2^). That is: 



(S(t 1 )S(t 2 )) = (S(t 1 -t 2 )) { -^l, (26) 

On 
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P\Hi t2) — i \*() 

V[0«(— 2i, *x|0) -0 R (-i,ti|O) 2 ] [fo(-2i,* a |0) - (j> R {-iMW\ 

3.2. The stock price differences 

The stock price differences are actually an usual variable involved in data analysis. We 
can define the stock price difference by 

Z{t-T) = S{t + T)-S{t) } (28) 

or the relative stock price difference in the following way 

Y(t;T)s s( t+ ^-s( t) (29) 

In consequence, we can do a similar analysis to these variable as the one already 
done for the stock. We want to obtain the several expressions for these variables in 
terms of the return characteristic function. 

For instance, the pdf of Z(t; r) is related to the stock price pdf as 

/■oo poo 

p z (z,t;r)= ds ds' S[z - (s - s')] p s (s, t + r; s', t\S ), 
Jo Jo 

and the characteristic function is 

(f> z (u,t',r) = ds d8'^ a - a ^ P s(8,t + T;8',t\So) 
Jo Jo 

/oo 
dr ,/ e iu,S (expr-expr ) pR ^ f + r; ^ £ |g) 
-oo 

/oo roo 
dr / dr' j»So(expr-exvr ) ^ _ ^ t |q) ^ 
-oo J — oo 

where Markovian and homogeneity conditions have been implemented. Now, as before, 
we can expand the exponentials in order to give an expression in terms of the return 
characteristic function. Thus, 



n=0 711 k=0\ K / J ^ ° 

/OO 
dre k( - r - r) p R (r-r',r\0) 
-oo 

= E til) (-V k M-in, t\0) M-ik, r|0). (30) 

t^o n - k=o \ k J 

Hence, its first moment is obtained with the definition in equation (0) 

(Z(t; r)> = S <M-i, t\0) r\0) - 1], (31) 

but the second moment needs a longer calculation with the following final expression 
(Z(t; r) 2 ) = S 2 0«(-2i, t\0) [1 + <f> R (-2i, r|0) - 2 <f> R (-i, r|0)]. (32) 
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Var[Z(t; r)] = S 2 {fo(-2i, *|0) [1 + <f> R (-2i, r|0) - 2 <f> R (-i, r|0)] 

-<M-i,t|0) 2 [l + i? (-i,r|O) 2 -2 i? (-i,r|O)]}. (33) 

Finally, we can study the correlation between stock differences and the stock itself. The 
correlation can be derived with equation ( f25|) and results 

(Z(t;r)S(t)) = (S(t + r)S(t)) - (S(t) 2 ) = S 2 <M-2M|0)[<M-i,r|0) - 1]. 

Simple manipulations that take into account equation ([H]) let us write 

(Z(t;r)S(t)} = ^^(Z(f,T)}. (34) 

And therefore the correlation coefficient is nonzero as it happens with the correlation 
coefficient between W(t; r) and R(t) presented in equation (|i"4|) . 

On the other hand, we can perform an equivalent analysis but for the relative stock 
difference, Y(t;r), given by equation (|29|) . Thus the characteristic function is 

M",t;T)=e- i " ^ t-L<f> R (-i n , r|0). (35) 

Observe that characteristic function for Y(t; r) is very similar to the one for the 
stock (pi)) but now the function is evaluated at time r, So does not appear, and an 
imaginary exponential appears in equation (j35D - This is not surprising since, due to the 
independency on time t in equation (^), the variable Y can be expressed as 

Y(t) = S(t)/S - 1 

thus implying the properties above mentioned. First and second moments, and variance 
are 

(r(t;r)) = fl (-i,r|O)-l, (36) 
(Y(t; r) 2 ) = 1 + <M-2i, t\0) - 2 ^(-i, r|0), (37) 
Var[F (t; r)] = fl (-2i, r|0) - ^(-i, r|0) 2 . (38) 

In addition, the correlation between Y(t; r) and S(t) is 

(Y(f,T)S(t)) = (Z(f,T)), (39) 

but from equations (|3lD , and ( ^BD we get 

(Z(t;r)) = (S(t))(F(t;r)). (40) 

Therefore, S(t) and r) are uncorrelated stochastic variables. 
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We take into account the results of the two last subsections and characteristic function 
of the return for the Wiener case of equation 

For the stock itself we can derive the following expressions. Therefore, we will 
subsequently obtain the first moment for the stock 

(S(t)) = S e^^ 2 ) 4 , (41) 



where it can be easily seen from equations ( 41|) and (pi]) that 



111 



s 



o 



( fi + a 2 /2)tjt f it= (R(t)), 



due to what is called the spurious drift component. Indeed, we can obtain the second 
moment which is 

(S(t) 2 } = S 2 e 2 ^ +a2)t . (42) 
In addition, it is possible to express the characteristic function as follows 

M^t\So) = t ( -^^ 2 e^, (43) 



n=0 



nl 



and the correlation function for the stock as 

(SitJSfo)) = S 2 e(^+^ 2 )^- te ) e 2 ^+- 2 )*2 

= g2 e M*i+*2) e | CT 2 (ti+3i 2 ) 



(44) 



where we take into account equations (p^D and ([41])-(|4"2|). And, on the other hand, we 
have the characteristic function for the price differences 



y-^ { lu}l= >0) c (na) 2 t/2 e fint ^ 



n=0 



k=0 



n 
k 



'_l) fc e ( fc(J ) 2 */2 e /^_ (45) 



And we can derive the first moment 

(Z(t; r)> = S e^ 2 / 2 )* [ e ^ 2 / 2 ^ - 1], 
and second moment 

(Z(t; r) 2 ) = So e 2 ^^ [1 + e 2{ ^ T - 2e^^ T }. 



(46) 



(47) 



We can also make similar calculus but for the relative difference. The characteristic 
function is almost the same as the one for the stock: 

oo 



4>y(uJ,t\S ) 



E- 

n=0 



1UJ 



111 



and the first and second moments are respectively 
(y(t;r)> =e( fx+a2/2 ) T - 1, 

(Y(t; t) 2 ) = 1 + e ^ +a2 ) 2r - 2 e^^ 2 ) 1 



(48) 

(49) 
(50) 
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The first moment of Y(t; r) is directly related to the first moment of the stock S(t) but 
differs from the W(t; r). However, when r is small, we thus have (up to first order) 

(Y(t; r)> ~ (// + a 2 /2) r^^r = (W(t; r)>, 

and, similarly, the second moment 

(Y(t;r) 2 } ~ (W^(t;r) 2 ) ~ <x 2 r. 

The main difference between the two averages is the spurious drift component which 
appears in the first moment of Y(t; r). We do not present a similar analysis on averages 
over S(t), since calculus leads us to complex expressions depending on time t. This 
feature also present in the Z(t;r) case and we will show in Section 5 that this fact is 
empirically observed in real markets. 



4. Empirical probability distributions 

Tick data from the Standard & Poor's 500 stock cash index is a good example of a 
multiplicative process, and shows the differences between working in terms of the return 
W(t;r) and in terms of the stock Z(t;r). In this case, time series analysis is able to 
give the empirical probability distributions of the differences between variables which 
are evaluated at two distinct times. 

For instance, stock price differences defined as Z(t; r) = S(t + r) — S(t) are plotted 
in Figure [I]. The two graphs give the probability distributions in tick data units when 
t = 1 minute. The first graph shows one-minute stock differences for the seventeen 
different years, from 1983 to 1999. And the second one only shows stock differences 
pdfs for the years: 1983, 1984, 1986, 1990, 1998. In this second plot we can see how the 
wings become fatter as the time increment t, in years, increases exponentially. 

Data analysis with time series assumes an annual periodicity, in the sense that 
Z(t; r) ~ Z(t + A; r) (in probability) where A = 1 year. In fact, we plot Z ignoring the 
time t variation and fixing r equals to 1 minute. As we have seen from equation (|30|) , 
the probability distribution depends on time t, fact which is in contradiction with the 
assumption that Z(t; r) is (in statistical sense) similar to Z(t + A; r). We see in figure [l] 
that plots become fatter as time t in years increases. 

On the other hand, Figure |] shows the return differences defined with the function 
W(t\ t) = R{t + r) — R(t). We plot the same type of probability distributions of the 
stock differences but for the return differences. If we compare the two figures we see 
that, in the returns case, the behavior does not change dramatically over different years 
(i.e., with time t) as the case of the stock differences. In effect, as it is proved in 



equation fll3]) , the return difference, when pdf is that of equation fll~5l), does not depend 
on time t and W(t; r) = W(r). 

We could also have plotted the relative stock differences but we do not think to 
be necessary. Relative difference Y is a better estimating variable than Z but it is not 
stationary unless we modify their value in an addequate way. We will study this in the 
next section. 
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-100 -50 50 100 

Stock differences (1 min.) 




Figure 1. We show the empirical pdf's for the tick data of the Standard & Poor's 500 
stock cash index differences. First graph involves one-minute stock differences for years 
ranging from 1983 to 1999. Second graph is a detail of the previous graph plotting 
pdf's of for years exponentially distributed between 1983 and 1998. 
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Return differences (1 min.) 



Figure 2. We show the empirical pdf's for the return differences of the Standard 
& Poor's 500 cash index. First plot involves one-minute returns differences for years 
ranging from 1983 to 1999. Second graph is a detail of the previous graph plotting 
pdf's of years exponentially distributed between 1983 and 1998. 
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5. Data analysis and estimators 

This section goes deeper in the study and comparison of W(t; r) and Z(t; r) estimators. 
Section 2 and 3 give us all necessary tools for measuring their quality, and thus giving 
correctly the first and second moment of the stock and return stochastic variables. In 
this way, we will also study some facts of Y(t; r) estimator and see that in some sense 
this estimator is halfway between the stock and return differences. 



5.1. Estimators for the first moment of the return and stock 

We define the following two sums for estimating the first moment: 

1 N ~ x 1 
M w (t; r,N) = - £ W(t + nr; r) = - W(t; Nt), (51) 

iV n=0 iV 
1 N ~ X 1 

M z (t; r, N) = - £ Z(t + nr; r) = - Z(t; Nr). (52) 

iV n=0 iV 

In latter expressions, we have synthesized the sum with the definition of the stock 
differences given respectively in equations (12) and (|2"8|). We want to study the quality 
of those estimators 0. Their averages are 

(M z (t;r,N)) = ^ (Z(t;Nr)), and (M w (t;r,N)) = ^ (W(t;Nr)). (53) 

We can easily derive the average over M w if we take into account equation (Tl3|). 
We see there that W evolves in the same way as the return and thus 

(M w (t;r,N)) = ^ (R(Nr)). 

However, this first moment can be decomposed in a sum of N equivalent terms although 
evaluated at several different times — see equation (ffD). That is: 

(M w (t-r,N)) = (R(r)). (54) 

Hence, the average of Mw estimator is equal to the first moment for the return. In the 
first plot of figure ^|, we have Myy estimator for the case when Nr = 1 year and r = 1 
minute. We observe that first moment changes from one year to another but with any 
specific trend. 

On the other hand, the average over Mz is obtained according to the 
expression (|3lj) , and reads 

(M z (fr,N)) = Is <M-M|o)[<M-i,iVT|o) - 1]. 

And, from equation (pl|), we finally have 



(Mz(t;r,N)) = { -§§ll(S(Nr)) - S ]. (55) 

In fact, the estimator is supposed to approach to the following value as the number of 
sample data increases 

lim (M z (t;r,N)} — > ^ ((S(r)} - S ), (56) 
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Figure 3. First and second moments of the one-minute return differences. We plot 
respectively M w (t; r = 1,T = 1) and [V w (t;T = 1,T = l)] 1 / 2 as a function of time t 
in years, from 1983 to 1999. Those functions are defined in equations ( pi] ) and (|67|). 
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Figure 4. First moment estimation with the one-minute stock difference. We plot 
Z(t;Nr = 1 year) /TV, for estimating (Z(t;r = 1 day)), as a function of time t (in 
years) from 1983 to 1999. Exponential growth with time t is observed as is also shown 
in equation (|60|). 



which when r is small (keeping only first order contribution) becomes 

lim (M z (t; r, AO) — (S(t)) (R(r)). (57) 

However, note that equations ( |5Tj| ) and ( |5TD will be valid only in case that following and 
equivalent limits are true 

^ ( i,jV T |0)-l <SW>-S, 
AT-+0O TV Lt|0) — 1] AT->oo j\T [(£(t)} — So] v 7 

Unfortunately this is not true in general. For instance, the Wiener process has the 
following limits 

e (A»+cr 2 /2)Nt _ i 

lim — —. ■ 9 , n , — > oo, (59) 

where we take into account equation (^TJ). Although r is very small, the term will tend 
to infinity as A^ approaches to infinity. Therefore, for this case, the stock differences 
estimator of the first moment is a biassed and not consistent estimator [0. 

We also observe the (Mz) depends on time t and, for the Wiener case, estimator 
evolves in average as 

(M z {t; r, AO) = S e^ 2 ^ -j- [ e ^ 2 ^ - l] , (60) 

which is derived taking into account equations (flop and (^5|). We see that the average 
grows exponentially with time t. This phenomena is also empirically observed in figure f| 
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for the one-minute stock differences graph. In this case, r = 1 minute, Nt = T = 1 
year and t = kT is evaluated in years. Therefore, equation ( |60"D reads 

(M z (k;r = lmin., N)) = S e ^+° 2 / 2 ) kT 1 



>W/2)T _ y 

where we can see that average grows exponentially with k similarly to figure f| which 
plots the first moment estimator in terms of k from 1988 to 1999. 

We may now study the limiting value of N for which Z gives a good estimation of 
(S(t)). For this to be possible 

N[(S(t))-S ]~ ' 
which is equivalent to demand 

(Z(t;NT))~N(Z(t,T)). (61) 

Assuming that the market follows the Wiener process, it can be shown that M z gives a 
good estimation for (S(t)) when 



-(/i + a 2 /2)iVr«l. 



(62) 



In addition, the variance of the estimator will determine us whether is an efficient 
estimator or not. Thus, 

1 



Vax[M z (t;T,N)] = — 2 [(Z 2 (t;Nr)) - (Z(t,Nr)Y 



(63) 

Taking into account equations (j3~I|) and (|32|), we finally obtain the variance of the 

estimator in terms of the characteristic function 
1 



C2 
O 



Var[M z (t; r, N)] = <j> R (-2i, t\0) [1 + fl (-2i, Nt\0) - 2 Nr\0)) 



-0|(-i,t|O) [<j) R (—i, Nt\0) — l] 2 . 
This can be represented in terms of the moments of the stock as 



(64) 



Var[M z (t; r,N)} 



A^ 2 Sn 2 



S 2 + (5(iVr) 2 ) - 2 S (S(iVr)) 



[(5(iVr))-5 ] 2 . 



And it is said to be a good (i.e., efficient) estimator when its variance tends to zero as 
A" tends to infinity [0. 

For the particular case of the Wiener process we will have 



V aT [M z (t;r,N)} = J ^{ 



1 + e 2(n+* 2 )Nr _ 2 e (n+a 2 /2)Nr 



^+a 2 /2)t 



1 + e 2( M +<r 2 /2)Nt _ 2 e 2{^+u 2 /2)NT Y 



(65) 



which also diverges as A" — > oo. We then conclude that, in general, this estimator is not 
efficient. For the Wiener case, the estimator is efficient only in case that A" is limited 
by the maximum value given by equation fl62f ). 
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Figure 5. Square root for the variance of the one-minute stock differences. We here 
plot [Vz(t;r — 1, TV)] 1 / 2 as a function of time t in years, ranging from 1983 to 1999. 
Exponential growth with time t is observed as is also shown in equation wW. 



5.2. The variances for the stock and the return differences 
We define the estimator for the variance Z(t; r) and W(t; r) by 

i N-l 



V z (t;r,N) 



V w {t;r,N) 



N 



- ^[Z(t + nr;r)-M z (t,r;iV)] 2 . 



n=0 
N-l 



N 



^2[W(t + nr;r)-M w (t,T;N)} 2 . 

J- „ — n 



(66) 



(67) 



n=0 



Similarly to the case above, we average the variance estimator in these two cases. Thus 
we have 

N-l _ i 



1 IM—l 

(Vw(t; r, AO) = £ [Wt + nr, r) 2 ) 



and 



n=0 



iV-1 



iV(AT- 1) 



{Wit-NrY), 



(Vz(t;r,N)) 



N 



- £ [(Z(t + nr,r) 2 ) 

1 n=0 



N(N — 1) 



(Z(t;iVr) 2 ), 



where we have used equations (|5T| ) and (|52|) . For the return differences case defined in 
equation (|12| ) and using the properties summarized in equation (|T3|), we obtain 



(W(t;r,iV)) 



AT — 1 



(68) 
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But if we particularize to the Wiener case, the variance average is directly related to 
the variance of the return, i.e., 



a 2 r. 



(V w (t;r,N)) 

where we take into account equation (|H|). 

And for the stock differences, we find similar divergences and limiting values for 
sample data to that of the case of the first moment estimator M z . Let us show this. 
We can write in terms of the characteristic function of the return with the help of 
equation (]32[) the average variance defined above 



(V z (t;r,N)) 



C2 
°0 



N- 



N-l 

+ <M-2i, r|0) - 2<M-i, r|0)] £ fo(-2i, t + nr\0) 

1 ' n=0 



i[l + fo(-2i, Ar|0) - 2fo(-i, Mr|O)]0 fl (-2i, t|0) 



Moreover, in terms of averages over the stock given in equations (^l|) and (|22j) , we have 
that the average is 

1 



(V z (t;r,N)) 



N-l 



[s 2 + (s(Tr}-2s (s(T)}} N ^^ S{f : nryl) 

n=0 



02 
O 



1 



[So 2 + (^(^)) 2 - 25 (S(Nr))) 



(S(tf) 



For the case when market model is the Wiener process, we then have 



(V z (t;r,N)) 



N-l 



r 2(LL+v 2 )t 



1 



e 2(n+cr 2 )Nr _ ^ 



N 



e 2( At +a2)r _ I 
1 _|_ e 2( M + ( r 2 )Vr _ 2e (M+<r 2 /2)iVT 



(69) 



Analogously to the first moment A tends to infinity the average of the estimator 

diverges. And, the estimator will be valid only when N obeys condition (p2;). If this 
condition holds and keeping r small, equation (|69f) proves that V z is a good estimator 
for the volatility. Hence, 

(V z (t;r,N)) = S 2 e 2 ^ + ^ t a\ (70) 

where we also observe that this estimator also grows exponentially with time t. This 
phenomena is also empirically observed in figure [|. Observe that last expression can be 
rewritten in terms of the second moment for the stock given by equation ([42]), that is: 

(V z (t;r,N)) = (S 2 (t)) a 2 r. 

Hence, for avoiding this divergence with t, a possible solution is to consider the estimator 
in the following way 

(Vz(t;r,N)) 



(S 2 (t)) 



a 2 r. 
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6. Conclusions 

Data analysis in financial markets is a very important issue due to the strong demand 
of higher precision in the estimation of parameters describing markets dynamics. For 
this reason, we have studied the stock and return differences when the return process 
is Markovian and homogeneous. Starting from the return characteristic function, we 
have derived the first and second moments, the variance and the correlation for the 
return, the stock price, and the return and stock differences. We have also obtained 
these expressions for the particular case when prices are driven by a Wiener process. 
After these calculations, we have compared the data analysis performed with the return 
and stock prices differences in the particular case that data source is the Standard & 
Poor's 500 cash index. 

We have intended to stress the importance in the way we manage financial database. 
It is well-known that stock data follows a multiplicative stochastic process but in some 
situations was, and still is, preferred to handle stock differences instead of taking return 
differences, that is: S(t + r) — S(t) instead of R(t + r) —R(i) [p], ||. The usual reason for 
doing this is that when r is small one can approximate the logarithm differences with 
the stock differences. We have showed that in general it is not true since estimators for 
the stock differences are biassed and not efficient. The approximation will be valid only 
when the sample data is smaller than a certain "critical" value and we have obtained a 
rule for the estimation of this value. 
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